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The theoretical analysis and numerical simulations of the magnetron operation with 
a feedback loop were performed assuming that the delay of the electromagnetic wave 
propagating in the loop is constant whereas the phase of the complex feedback reflec- 
tion coefficient is varied. Results of simulations showed that by a proper adjustment 
of values of the time delay and phase of reflection coefficient that determines phase 
matching between the waves in the resonator and feedback loop, one can increase the 
magnetron's output power significantly without any other additional measures. 
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I. INTRODUCTION 

Magnetron oscillators are widely used sources of microwave radiation operating over a 
wide range of powers and frequencies and having various applications!. A typical feature 
of the magnetron operation is a relatively broad and noisy spectrum 2 . However, for many 
applications (radar, communication systems, charged particles accelerators, etc.), it is very 
necessary both to narrow the frequency spectrum and to reduce the noise. An improvements 
in the microwave spectrum can be achieved either by using a special magnetron design, which 
makes starting conditions more preferable for the operating mode, or by feeding external 
monochromatic electromagnetic radiation into the magnetron cavity. In the latter case, 
this controlling radiation stabilizes the magnetron operation similarly to other nonlinear 
oscillators-. In fact, the injected electromagnetic radiation can be produced by either an 
external source [injection-locked magnetron (ILM)£~— ] or by the same magnetron [self-inject- 
locked mag netron (SILM) 8 ]. In the latter part of the magnetron output power is 

injected back into the magnetron cavity via a feedback loop. It is worth noting that this 
feedback loop can appear as a result of the magnetron radiation reflection from the load. It 
had been shown« that ILMs and SILMs are able to fix the operating frequency and phase, 
and essentially decrease the noise of the radiation. 

Recent investigations*^ show that the efficiency of the relativistic S-band magnetron 
can be significantly (up to 40%) increased when part of the radiated electromagnetic power 
is reflected into the magnetron cavity. This effect, to our knowledge, has not yet been 
discussed. In this paper we suggest an explanation for this effect using a model of the 
magnetron which considers it as an active nonlinear ring resonator. This model can be 
applied also to coupled ring resonator-fiber systems which are of interest for optoelectronics 
and communication (see, e.g.,— >^). 

II. GEOMETRY OF THE PROBLEM 

In the model, a magnetron is considered as a ring resonator (see Fig. [T])with double 
degenerated eigenmodes, clockwise |+) and anti-clockwise |— ), propagating waves having 
equal frequencies uj. Only one wave (mode |+), for definiteness) can be excited due to 
the resonant interaction with the rotating electron flow when the electron azimuthal drift 
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velocity is close to the wave phase velocity. The energy stored in the |+) wave (hereafter 
active mode) leaks partially from the resonator to the waveguide (see Fig. [T^i). 

The output end of the waveguide is connected to a certain payload, for instance, an 
antenna. It is assumed that the matching between the waveguide and the load is not ideal, 
i.e., part of the wave energy is reflected and enters back into the magnetron cavity. This 
reflected wave excites two contra-propagating waves, |+) and |— ), in the magnetron cavity. 
The active, |+), wave is amplified by the electron flow, whereas the passive wave, | — ), 
propagates without amplification. Thus, one can neglect the passive wave and consider 
the waves' paths as shown in Fig. [Tb. By convention, the paths of rightward and leftward 
propagating waveguide modes (| — >) and | <— ) modes) can be considered as separated paths. 
Namely, the wave which is reflected from the waveguide output end comes back to the 
resonator along a feedback loop (see Fig. [lb). In this scheme, the ring resonator is side- 
coupled to the waveguide. The coupling between the resonator and waveguide modes is 
similar to a phase matched coupling. This means that only the co-propagating waves are 
coupled: modes | — >) and |+) (modes | <—) and |— )) are coupled, whereas modes | — >■) and 
|— ) (modes | <—) and |+)) are not coupled. This configuration, i.e., side-coupled waveguide 
and ring (or disk) resonator, has widely been explored in photonics (see, e.g.,— ^). The 
general properties of this system were reported in Ref.— and the results of this work will be 
used in the next section. 

III. SMALL SIGNAL MODEL 

Let us consider a supplementary problem, namely wave propagation in the waveguide 
coupled with a ring resonator (see Fig. |2]). Let ipi n ^ be the complex amplitude of the 
monochromatic rightward propagating waveguide mode ipin{t) — V'm^e - ^* on the left side 
with respect to the waveguide-resonator coupling cross-section. The complex amplitude 
Vw,u; °f the transmitted wave ipout(t) = r ^cmt,^~ wt on the right side with respect to this 
cross-section and amplitude are related as: 



*(w)V>»n,w, 



(1) 



where the transmission coefficient t(u) is defined as^: 



t(u) = 1 - i 



1 



L 
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(2) 



uj — uj r + iT Vg 



r,w 
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FIG. 1. Magnetron (ring resonator) and waveguide, (a) - Clockwise |+) (yellow bold arrow) and 
anti-clockwise |— )) (red bold arrow) propagating eigenmodes. Direct (dashed line) and reflected 
(solid line) waves in the waveguide, (b) - Ring resonator side-coupled to a waveguide with a 
feedback loop. 




FIG. 2. Auxiliary problem. Phase matched coupling between waveguide and waveguide modes. 
Notation are the same as in Fig. Q] 

Here u r is the eigenfrequency of the ring resonator, r = + r c , the coefficient is 
related to the resonator dissipative losses, the coefficient T c = \ V r ^ w \ 2 L/2v g is related to the 
energy leakage from the resonator because of coupling with a waveguide characterized by 
the coupling coefficient V r>w , L is the resonator perimeter, v g is the group velocity of the 
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resonator eigenmode. Eq. (j2]) is valid when the coupling is weak: 

V c /U r < 1. (3) 

Parameters T, Yd, and T c can be written in terms of the resonator Q-factors: F = 
u r /2Qtot, T d = u r /2Q diss , and r c = u r /2Q leak . Here Q tot is the total Q-factor, and Q diss 
and Qieak are the Q-factors responsible for the energy dissipation and leakage, respectively. A 
negative dissipative Q-factor, Qdiss < 0, corresponds to an active resonator, which amplifies 
an incoming wave. 

The output wave amplitude Vw,w is a superposition of two amplitudes, namely, the 
amplitude of the wave which crosses the coupling region directly through the waveguide, 
and the amplitude of the wave that enters the waveguide from the resonator. The two terms 
in the right side of Eq. §2§ describe these two paths of the wave transmission. The first term, 
unity, describes the wave's direct transmission through the waveguide^; the second term 
describes the wave's transmission through the resonator. Thus, the output wave amplitude 
ipout,ui is the sum of two amplitudes: 

i>out,u = t(L0)lp iTlt u = 1pdirect,u> + *Pres,uj, (4) 

where ipdirect,^ — 4>in,uj is the amplitude of the wave transmitted directly through the waveg- 
uide, and ipres,Lj is the amplitude of the wave that transmitted through the resonator: 

Aes,u> = [t(u) - = ; -r^ x - VW>- ( 5 ) 

Eq. ()5]) is the Fourier transformation of the differential equation: 

d ^res(t) , . U r f 1 1 \ UJ r 

— + tU r 1p res {t) + — — + — Aes{t) = -pr Winit)- (6) 

Thus, one can write the following equation for input and output waves having arbitrary 
shapes: 

tpout(t) = tpin{t) + Aes{t), (7) 

where ip res {t) is the solution of Eq. ([H]). Here, let us remind that the wave ipres{t) is the wave 
that leaks from the resonator into the waveguide. The wave field in the resonator can be 
significantly larger than the wave field ip res (t) in the waveguide [weak coupling between the 
resonator and waveguide, see Eq. However, these wave fields are linearly proportional 
to each other and the field ip r es{t) can be considered as a normalized field in the resonator. 



A supplementary problem that is being considered helps one to describe the system with 
the presence of a feedback loop (see Fig. [lb). In a general case, the wave propagation along 
the feedback loop is described by a functional equation: 

1) in (t) = F{^ out (t)}, (8) 

where F is a certain functional. In the simplest case of a non-dispersive feedback loop this 
relation has the form: 

ip in (t) = qip out (t-T), (9) 

where T is the time of the wave propagation along the loop (delay time), and q is the complex 
feedback coefficient, \q\ < 1. Below we will consider this simple model only. 

Eqs. (|BD, (J7J), and © [dED in the general case] being supplemented by initial conditions, 
are a closed system that describes the small signal theory of a magnetron with a delayed feed- 
back. Asymptotic solutions of these equations are rather simple. There is an exponentially 
either growing or decaying output wave amplitude, depending on whether the system's self- 
excitation threshold is exceeded or not. In order to describe the wave amplitude saturation, 
non-linear terms should be incorporated in these equations. 

IV. NONLINEAR MODEL 

Without the feedback loop, i.e., when q = (ip in = 0), Eq. (j6j describes a linear resonator 
with the self-excitation threshold defined as: 

Qdiss + Qleak = 0. (10) 

When Qdl 3S + Qi^k < 0, the field amplitude ifj re8 (t) growths exponentially. The increase 
in the amplitude can be limited by the replacement of the linear dissipation term Q^ iss < 
in Eq. © by the nonlinear one: Q^ ss ->■ Q^ ss = Q^L i 1 ~ \ Aes(t) | 2 hi), where 7 d is the 
parameter that is related to the saturation amplitude. With nonlinear Q~^ s Eq. (jBJ takes 
the form: 

VVes(t) = 0. (11) 

This is the Van der Pol equation which is widely used as a model of a nonlinear auto- 
oscillating dynamical system^. Another effect which, in general, should be taken into ac- 
count, is a nonlinear shift of the resonator eigenfrequency. This effect can be incorporated 
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for in the model by the replacement in Eq. fl6]) u r — > u r (1 — \ip res {t)\ 2 / ^ f) , where 7/ is the 
parameter that is related to the frequency shift. The latter results in describing the Duffing 
oscillator. Both Van der Pol and Duffing models are typical examples of nonlinear two- 
dimensional oscillating systems^. Duffng's model describes non-isochronal oscillations, and 
Van der Pol's model describes a limit cycle when the equilibrium point loses its stability. It 
is reasonable to suppose that, due to their generality, these models can be used to describe 
the properties of the system under consideration. To summarize, let us present below the 
equations which will be analyzed: 

dijj res (t) 



dt 



+ iu r (l - \Aes(t)\ 2 /lf) Aes(t) + 



0J r 



Q 



1 1 
— + 



leak 
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l-\lP res (t)\ 2 /% 



-4>in(t), 



1pres(t) = 

^Cleak 

il>in(t) = q^out(t - T), 

1pout(t) = ll>in(t) + Aesit). 



(12) 
(13) 
(14) 



V. ANALYSIS 

One can rewrite Eqs. f|T2|) - f|T4|) in the dimensionless form using dimensionless time r = u r t 



and normalized amplitudes Aj(r) = V>i(i~)e lT 7 d V Qleak/ (Qleak \Qdiss\)v 

- ^|QcKss| |^4res(' r )| A res (r) ~ 

CLT 7j 

-k [l - |A, es (r)| 2 ] Aesir) = --^—A in {r), (15) 

^ ^Cleak 

A in (r) = qe- ie A out (r-9), (16) 
A out {r) = A in {r) + A res {r), (17) 

Here 6 = u r T is the dimensionless delay time and k = {Qi ea k — I Q diss I) /Qieak\Qdiss\- 111 
Eq. ( fTBT) it was taken into account that Qdiss < 0. 

Eqs. ( fT5l) -( fT71) have monochromatic solutions with constant amplitudes. Setting Ai(r) = 
Bie~ iUT , where v = (oj — u r )/uj r is the normalized frequency shift, we arrive at the following 
set of algebraic equations: 

7^ 1 1 

%vB res ^ I Qdiss I |-Sres| B res — K (1 |-Bres| ) B res — - -Bjn; 
7/ 2 X ' Qleak 

B in = qe- ld+lv9 B ouU (19) 
B ou t — B in + B res . (20) 



Excluding B in , one can obtain: 

/ 2 \ i i -iO+ivO 

i [v + K\Q*,,A f f \B„,\^ + -« (i - = _ r A__. (21) 

Eq. ( l2Tj) determines both frequency shift v and normalized power in the resonator P = \B r ~ 
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(P = 1 without the feedback, q — 0) . In the simple case when the nonlinear dispersion absent 
(7/ — > 00) one can rewrite the complex equation ( )2T|) as: 
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(22) 



1 KQkafc 1 + |<?| 2 - 2|g| cos(6' - z/6» - a)' 

where a = argg the phase of the feedback coefficient. 

As it follows from Eq. ( 1231) . the normalized power can be varied in the range: 
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When the left boundary in (|24|) is negative, the minimal possible power is equal to zero. 

Near the self-excitation threshold, the parameter ^Qi ea k = (Qieak — \Qdiss\) /\Qdiss\ is 
small. In the latter case the feedback can noticeably increases the power P in the resonator. 
The inequalities (12^1) determine potentially-possible boundaries of the power variation. Upper 
and lower boundaries correspond to such values of the frequency shift v for which cos(6 — v9 — 
a) = =f1, respectively. However, at designated values 9 and a, the frequency v is determined 
by the solution of the dispersion equation f[2~2"l) . The proximity of the power value to the 
upper possible value depends on whether the solutions of the dispersion equation contain 
eigenfrequencies v for which cos(6 l — vQ — a) is rather close to —1. 

It is convenient to present the dispersion equation f l24l) as: 

0\q\ sin 09 . . 

V = e-a+ m * 25 

Qleak (1 + \q\ 2 ~ 2\q\ COS (p) 

where <p = 9 — a — Ou. If the feedback coefficient is small, \q\ 2 1, and 

e<Qiea k /\q\, (26) 

Eq. (1251) has one solution (p(a) (assuming delay time 9 = const). The inequality (|26|) ) will 
be considered as a short delay time condition. In the case of a reflection coefficient phase a 
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variation, the value cos[<^(a:)] varies in a full range from —1 to +1. Accordingly, the power 
P varies in a full range defined by Eq. (123]) . In the opposite case of a long delay time , when 

0^Qleak/\q\, (27) 

the dispersion equation ( 125]) has many roots, and the eigenfrequency v for which cos[# — 
v(p)Q — a] is close to —1 is always present. This means that Eqs. (|22|) and (|23|) always have 
solution for which the power P value is close to its upper limit defined by Eq. (EMI) . 

The analysis presented above shows that a variation in the reflection coefficient phase can 
produce a rather strong power variation, if the inequality (|26|) is satisfied. In the opposite 
case, when the inequality (1271) is satisfied, the power is close to its maximal possible value, 
which can noticeably exceed the power in the resonator without a feedback loop. 

VI. NUMERICAL SIMULATIONS 

Examples of numerical solutions of Eqs. (fT5T) - (fTTj) are presented in Figs. [3] and HI 




Normalized time 

FIG. 3. Dependence of the normalized power on time. "Short" delay time < Qi ea k/\q\- Red 
dashed line corresponds to q = 0. The curves are marked by corresponding a /it values. \Qdiss I = 15, 
Q leak = 20, q = 0.2, 9 = 70. 

Fig. [3] demonstrates how the power saturation level depends on the phase a when the 
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inequality ( 126|) is satisfied. Parameters used for numerical simulation were chosen close to 
their estimations for the experimental setup described in Refs.- 1 ^: |Q* SS | = 15, Qi ea k = 20, 
q = 0.2, 9 = 70. The results of simulations showed that the microwave generation in the 
magnetron can be either totally suppressed (a = 0) or enhanced by ~ 40% (a ~ 7r) as 
compared with the magnetron operation without feedback (q = 0). 
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FIG. 4. The same as in Fig. [3j "Long" feedback loop, 6 3> Qi ea k/\q\- \Qdiss\ = 15, Qi ea k = 20, 
q = 0.2, 9 = 200. 

Fig. H] is related to the "long" feedback loop when the inequality f T2T|) is satisfied. In 
agreement with the analysis presented above, the power saturation level is almost indepen- 
dent of the phase a and is close to its maximal value for the system with a "short" feedback 
loop (short delay time). 

VII. DISCUSSION 

The theoretical analysis and numerical simulations of the magnetron operation with a 
feedback loop were performed assuming that the time 9 of the wave propagating in the 
feedback loop is constant whereas the phase a of the reflection coefficient q is varied. These 
assumptions were used only to simplify the solution. In general, the solution of the dispersion 
equation fl25|) depends on both phase a and delay time 9. Variation in the value of 9 
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affects the solution in the same manner as the phase a variation. This similarity becomes 
important particularly when 9 ^> 7r and the condition fl26l) is satisfied. Thus, the length 
of the waveguide that connects the magnetron with a payload (an antenna, for instance) 
can affect the magnetron characteristics essentially. Results of simulations showed that by 
a proper adjustment of values of 9 and a that determines phase matching between the 
waves in the resonator and feedback loop, one can increase the magnetron's output power 
significantly without any other additional measures. 
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